Feedback effects and the self-consistent Thouless criterion 
of the attractive Hubbard model 
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We propose a fully microscopic theory of the anomalous normal state of the attractive Hubbard 
model in the low-density limit that accounts for propagator renormalization. Our analytical con- 
clusions, which focus on the thermodynamic instabilities contained in the self-consistent equations 
, associated with our formulation, have been verified by our comprehensive numerical study of the 

same equations. The resulting theory is found to contain no transitions at non-zero temperatures for 
all finite lattices, and we have confirmed, using our numerical studies, that this behaviour persists 
. in the thermodynamic limit for low-dimensional systems. 

i 



The anomalous normal state properties of the high temperature superconductors have, amongst other things, 
highlighted the limitations of present-day many-body physics. In spite of repeated attempts Jl|-[| to formulate a 
controlled fully microscopic many-body expansion, at present both microscopic and semi-phenomenological theories 
require vindication from Monte Carlo or some other "exact" procedure. 

To serve as further and concrete motivation for the work described below, consider the paper of Schmitt-Rink, 
Varma, and Ruckenstein Q|. They proposed, through a simple formulation of the problem, that in two dimensions 
O | (2d) any electronic effective two-particle interaction that was attractive would lead to a destabilization of the Fermi 
surface, and thus to non-Fermi liquid-like physics. One may consider this result as being based on the idea that 
' the Fermi surface is destroyed via the presence of two-particle bound states that "drain away" all electrons to the 
. bound state. The result is a superconducting instability which, at best, occurs only at zero temperature. The T = 
predictions of this model are complicated: they found a critical point, corresponding to all electronic densities, with all 
electrons paired — whether or not these paired electrons were in a phase coherent state was not addressed, although 
£Q ■ that is what one would expect. 

Many questions and criticisms of this approach have appeared [^). In particular, by including multiple scattering 
in the electron number equation, Serene |J claimed that the physics of Ref. Q was eviscerated, and that Fermi liquid 
behaviour was in fact robust. Actually, as we show below, at the multiple-scattering level of formulation the physics 
. of Ref. Q survives such an improvement. 

Another criticism of the result of Schmitt-Rink, Varma, and Ruckenstein Q is that Luttinger's theorem is not 
satisfied — this immediately suggests that renormalization of the single-particle propagators is important. As we show 
below, this is indeed the case, and, in fact, the physics of Ref. Q is dramatically changed when such renormalization, 
or self-consistency, is included. The problem with such a proposal is that the renormalization of these propagators is 
not a simple matter, and in this paper we address the question of the "appropriate" self-consistent theory, and apply 
our considerations to the attractive Hubbard model (AHM). 

The impact of self-consistency on the results of Ref. Q is related to the more general question of the study of feedback 
effects on the Thouless criterion. The Thouless criterion arises from the determination, with decreasing temperatures 
in the normal state, of the superconducting instability Such a formulation traditionally uses unrenormalized 
$_i ' electron propagators — see || for a particularly clear discussion of the physics that arises from such an approach. An 
important advantage of the Thouless approach is that it focuses on the non-superconducting state, so the possibility 
of non-Fermi liquid physics in the anomalous normal state can be studied simultaneously with the question of the 
superconducting instability. Thus, the issue of the self-consistent Thouless criterion can be addressed by repeating the 
analysis of, say, Ref. H, using renormalized propagators. This provides feedback, i.e. the single-electron propagators 
will "know" about the pairing tendency of the system as the temperature is lowered, and this in turn will influence, 
for better or for worse, that pairing tendency. 

What self-consistent formulation is appropriate? One suggestion comes from Levin and coworkers f9j|l0| , based on 
earlier work jlJ|T^-|T^] , in which only one specific combination of renormalized propagators (in the equations for both 
the pair susceptibility and the self energy) was proposed. Fresard, et al. (m), and then Haussmann ff5|| (and several 
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subsequent authors), examined the situation that arises when all propagators are renormalized. Another altogether 



different approach was adopted by Vilk and Tremblay 16 1, who used sum rules, symmetry, and other considerations, 



to construct a theory in 2d which required that renormalized propagators not be used. Our work is different from all 
of these attempts. 

For electron densities near half-filling there is no justification for omitting vertex corrections. However, the approach 
we adopt is designed specifically for low electron densities, for which we have demonstrated that vertex corrections 
are negligible. We find (i) no superconducting transition at non-zero temperatures for all finite lattices, curing a 
difficulty that is present with all mean-field theories; (ii) a self-consistent solution at T = that corresponds to a 
superconducting phase; and (iii) numerical solutions of the resulting self-consistent equations confirm our analytical 
results and demonstrate that in low dimensions our conclusions survive in the thermodynamic limit. Further, based 
on our numerical solution we predict the low-temperature form of the pair susceptibility. 

The above-mentioned papers examined a variety of different Hamiltonians. For concreteness, we used the attractive 
Hubbard model (AHM) on a d-dimensional hypercubic lattice, defined by 

H = -t c tn,<7 c m\cr - \U\ ^2n mA n md (1) 

(m,m'),o" m 

where the terms in this Hamiltonian have their usual meaning. 

When the electronic density is either close to zero ((n m ) ~ 0), or complete filling ((n m ) ss 2), the effective interaction 
is well described by the T- matrix approximation fl^ j . The T-matrix approximation corresponds to a self-energy given 
by (ignoring, for this publication, the Hartree term) 

U 2 

£(k,zw„) = -— Y 02(q,«fif) G a (-\s. + q,iQ,e-iuj n ) , (2) 

where G2 is the pair propagator, and G Q can be either Go (the non- interacting) or G (the renormalized) single-particle 
Green's function. The pair propagator is given by 

y 2 (q,«k) = z \~Fj~\ 7 — 7TT ' ( 3 

1 - \U\ x(q,i^) 

where the pair susceptibility, x, is a convolution of two single-particle propagators given by 

X(q, tilt) = G a > (k, iuj n )G all (q - k, iQ. t - iw n ) , (4) 

with G a i and G a » either Go or G. In these equations, (3 is the inverse temperature (ks = 1), iV is the number of 
lattice sites, and io n (Tie) are Fermi (Bose) Matsubara frequencies. Finally, Dyson's equation gives another relation 
between the self-energy and the renormalized Green's function, G. 

As outlined in Ref. Q , the Thouless criterion, signifying the instability of the normal state, corresponds to poles of 
the zero-centre-of-mass momentum retarded, real-time, pair propagator moving off the real axis into the upper half 
of the complex plane. From Eq. (|^) this simplifies to the condition (defining \ = \' + x") 

X '(q = 0,^ = 0) = l/\U\ , (5) 

(noting that x"(q = 0, ifle = 0) = 0). 

For future reference, let us make clear the instability that is found when propagator renormalization is ignored. 
Traditionally ||, one evaluates Eq. (|j) using bare propagators. The divergence of Eq. (3) then leads to an instability 
at a temperature T c which is identical to that predicted by BCS theory (e.g., see Ref. [8]). Below we discuss some of 
the effects of feedback on this scenario. 

As further information on the above no-feedback theory, note that researchers have begun to use the appearance of 
the pseudogap as an important test for the presence of strong correlations in model systems, and this requires that one 
obtain the dynamics, and not the static thermodynamic quantities, for the proposed model. However, to obtain the 
dynamics from even the simple non-self-consistent T-matrix theory involves some hard-to-control extrapolations to 
the real frequency axis. Recently, we formulated a partial fraction decomposition method, that is described elsewhere 
jUI , and which can be evaluated with a precision (with respect to pole locations and residues of the self energy) to 
a relative accuracy of 10 -80 . The resulting density of states shows a well-defined pseudogap at low temperatures 
pq|, very much reminiscent of experiment pjl. Further, our results for the single electron spectral functions derived 
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from the non-self-consistent T matrix are remarkably similar in character to the quantum Monte Carlo/ maximum 
entropy method spectral functions and density of states determined by the Sherbrooke collaboration [£0| (although 
their calculations are for a much higher electron filling). Thus, from a comparison to either experiment or "exact" 
Monte Carlo data, it appears that a T-matrix theory with unrenormalized propagators is a reasonable formalism with 
which to examine this problem. However, as discussed at the start of this paper (with regards to Ref . [Q ) , and as we 
now delineate below, this agreement is fortuitous — the predictions for the (static) thermodynamics for such a theory 
are in error, and thus the dynamics produced by such a formulation are of no value. 

To see how the thermodynamics predicted by Refs. Q^J make clear the failings of such a theory, consider the 
following: If one chooses a chemical potential that lies in the band (|/x| < W/2, where W is the single-particle 
non-interacting bandwidth), and then lowers the temperature, using the original Schmitt-Rink et al. (J] formulation 
applied to the 2d AHM, we have shown that the electron density actually diverges as the instability is approached 
(for constant a chemical potential). Use of the correct number equation, as follows from the number-conserving 
theory of Serene ||, leads to a number density of one at the instability temperature (see below). Then, this leads 
to the unphysical result that for any fixed electron density less than one, the chemical potential is always driven 
to half the bound state value (viz. below the band), indicative of the elimination/suppression of the Fermi surface. 
Further problems abound — while the electron density approaches unity at the instability, the two-particle correlation 
function (n m ^n m ^i) diverges — that is, this formulation has difficulties at "the two-particle level" This is clearly 
incorrect, since for the AHM it should be bounded between (n m ) 2 /4 and (n m )/2. 

We now demonstrate the above-mentioned behaviour analytically: We note that Eq. (|J) indicates that near the 
Thouless curve (T c ) the self-energy (Eq. (0)) is dominated by the q = 0, Qg — term. Thus for all finite lattices one 
finds 

vn, • ^ U2 X(0,0) r ( v ■ ^ 

£(k,«w n ) ~ - — — — — — — G (-k, -iw n ) . (6) 



N/31-\U\ x(0,0) 



Note the electron density equation is 
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(ftm) = 1 + J^nKe ( iuJ n ~ 6c - S(k, iUn))" 1 , (7) 

where £k is the band energy relative to the chemical potential, and that at the Thouless temperature the self-energy 
diverges (Eq. |J), and thus (n m ) = 1, as stated above Several example flow lines are shown in a (//, T) plot for 
a 2d square lattice in Fig. , in which it is seen that they approach the two-electron (hole) bound state energy below 
(above) the band as T — » 0. This results in a system of bound electrons (holes), and the Fermi surface (and hence 
the normal state Fermi liquid) is destroyed. 

To remedy these deficiencies within a fully microscopic diagrammatic theory, renormalized propagators are required. 
Unfortunately, there are a variety of different self-consistent T-matrix theories. These different theories arc based on 
the choice of a non-interacting Go or a fully renormalized G in the G Q 's that appear in Eqs. Based on the 

consideration of Marcelja |L2| and Patton jl3]], Levin and coworkers Hfio) have argued that the G a in Eq. (^) must 
be the non-interacting Green's function (to avoid states in the superconducting gap, which are absent in the weak 
coupling limit). While we agree with this "non-derivable" approach, we disagree with the claim that the appropriate 
combination of Green's functions in the pair susceptibility should be asymmetric (x ~ G Go). Firstly, these theories 
do not produce the full two-particle Green's function, but rather just the pair propagator. One cannot say that a 
pair propagator is conserving, but one could ask: Is the two-particle Green's function from which the pair propagator 
is formed conserving? It follows from Ref. (2| that a symmetric convolution (~ G G) leads to two-particle quantities 
that are conserving, and thus it is not required, nor necessarily wise, to use the asymmetric form (~ G Go) HJrJ for 
the pair propagator. 
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FIG. 1. Lines of constant density are plotted in the fi-T plane (each quantity shown in units of the bandwidth W) for a 
14 x 14 lattice. The contours shown here (from bottom to top) correspond to electron densities (n) = 0.1, 0.2, 0.3, . . . , 1.8, 1.9 
and are calculated for \U\/t = 4. The Thouless criterion line forms the boundary to the shaded superconducting region. Notice 
that all the density contours, except the (n) — 1 line, are deflected away from this area, and never cross into this area, as they 
proceed to lower temperatures. 



Based on these (and other |£2[) considerations, we examine the self-consistent theory given by Dyson's equation 
with G a = G and G a > = G a " = G in Eqs. (|,|). Following our earlier arguments, for non-zero temperatures and a 
finite lattice we can write the ansatz |23| ] near the superconducting instability 

A 2 

- J^W, ' (8) 

A 2 = U2 * ( °' 0) (9) 
Np(l-\U\ X (0,0)) ■ 1 1 

From Eq. (||) it follows that a T c > instability is indicated by A — > oo. However, evaluating x(0, 0) using Eq. (||) in 
Eq. (|) leads to 

*(°.°> - w 3 g • « 

and shows in fact, that x(0) 0) goes to zero in the limit of A — > oo. Thus, there is no self-consistent solution at non- 
zero temperatures for any finite lattice — this cures an obvious problem that exists for mean-field (non-self-consistent 
T- matrix) theories p4|. 

The only self-consistent solution of these equations occurs at T = 0. That is, by solving 1 = |£7|x(0, 0) in the T — > 
limit we obtain a modified gap equation given by 

i = H ( * ; 1 ^ ) (n) 

One finds that in the small \U\/t limit, A ~ exp(— 4W/|J7|); this shows, in contrast to what was claimed in Ref. [p|,^0[, 
that when both propagators in the pair susceptibility are renormalized one still recovers the weak coupling BCS limit 
p5[. Indeed, the form of the density of states at T = 0, 
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aah~aa°(o)^==0(^- a2 ) . (i 2 ) 

V uj z — A z 

is that of weak coupling BCS theory, and shows a fully developed gap. 

In order to further analyze the predictions of our theory, we have conducted a comprehensive numerical study 
of the resulting self-consistent equations, and have solved these equations down to the lowest temperatures that 
we can access (typically, these temperatures are l/100th of the mean-field BCS temperature). Further, we have 
obtained convergence of all static quantities to better than 1 part in 10 5 (at low temperatures this requires that our 
Matsubara frequency sums be performed with a frequency cutoff of AW, and we have thus used this same cutoff at 
all temperatures). An example set of data is shown in Figure 2 for d = 1; lattices from 8 x 1 to 128 x 1 were studied, 
and it was found that for L x 1 lattices, for L > 64 no changes with lattice size were encountered. Thus, these data 
represent the quantitative predictions of our renormalized propagator theory for one dimension in the thermodynamic 
limit (with the shown parameter set). As is clear from this Figure, the T — > limit of x(0, 0) is 1/\U\ (in units of 
1/t); however, x(0, 0) remains less than 1/|C/|, and thus no instabilities are encountered, for all nonzero temperatures. 
Clearly, as is predicted by our analytic theory, there is only a zero temperature transition to a superconducting state. 




1 ' ' ' 1 
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FIG. 2. x(T) = x(q = 0,i£le = 0) vs. temperature, T, both in units of 1/t, for the parameters indicated. The Thouless 
temperature, T-rh, in the (unrenormalized) non-self-consistent T-matrix theory (Go -Go) is indicated with an arrow; by including 
renormalization, this temperature is driven to zero (curve indicated with G • G). We have plotted results for 32 x 1, 64 x 1 and 
128 x 1 lattices; the results are essentially indistinguishable (and are in fact plotted on top of one another), and all extrapolate 
to 1/\U\ — 0.5/t at T/t — 0, to within an accuracy of better than 0.01%. 



Similar data and similar trends are found in two dimensions. However, if we impose the same strict convergence 
criterion that we did for our Id results (for systems in the thermodynamic limit) we find that we cannot faithfully 
extrapolate to the thermodynamic limit in 2d. For example, using a constant number of Matsubara frequencies leads 
to converged solutions which are unstable with respect to the inclusion of more frequencies, while using a fixed energy 
range (as we did for Id) leads to limitations in our computing power. (At present we are exploring different strategies 
to overcome this restriction.) Thus, the behaviour of this formulation in the thermodynamic limit in two dimensions 
is an open problem, and whether or not these numerics confirm the interesting Kosterlitz-Thouless prediction of Ref. 
pG| is unknown. 

In conclusion, we have formulated a self-consistent theory with renormalized propagators of the AHM. We employ a 
symmetric pair propopagator formed with renormalized single-particle Green's functions G, but "close" the self-energy 
with an unrenormalized Go- We have found a superconducting solution at T = 0, but for all finite lattices we have 
proven analytically and confirmed numerically that no non-zero temperature instabilities exist. The resulting zero- 
temperature superconductor has a conventional density of states in the weak coupling limit. Although it is unclear if 
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such a theory accounts for critical effects at T = 0, it is clear that this formulation represents a very useful starting 
point of a fully microscopic theory for the evaluation of, e.g., the spectral functions at non-zero temperatures (albeit 
for low electronic densities). 
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